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O ■ ABSTRACT 

O^l , We present parameter estimation forecasts for present and future 3D cosmic shear 

^ ■ surveys. We demonstrate in particular that, in conjunction with results from cosmic 

Q ' microwave background (CMB) experiments, the properties of dark energy can be es- 

timated with very high precision with large-scale, fully 3D weak lensing surveys. In 
particular, a 5-band, 10,000 square degree ground-based survey of galaxies to a me- 
dian redshift of Zm = 0.7 could achieve l-tr marginal statistical errors, in combination 
with the constraints expected from the CMB Planck Surveyor, of Awq = 0.108 and 
AiUa = 0.099. We parameterize the redshift evolution of w by w{a) — + Wa(l — a) 
where a is the scale factor. Such a survey is achievable with a wide-field camera on a 4 
metre class telescope. The error on the value of w at an intermediate pivot redshift of 
' z = 0.368 is constrained to Aw(z = 0.368) = 0.0175. We compare and combine the 3D 

, weak lensing constraints with the cosmological and dark energy parameters measured 

\^ • from planned Baryon Acoustic Oscillation (BAO) and supernova Type la experiments, 

' and find that 3D weak lensing significantly improves the marginalized errors on wq and 

. Wa in combination, and provides constraints on at a unique redshift through the 

' lensing effect. A combination of 3D weak lensing, CMB and BAO experiments could 

\ achieve Awq = 0.037 and Awa = 0.099. We also show how our results can be scaled 

Q^. to other telescopes and survey designs. Fully 3D weak shear analysis avoids the loss 

I ' of information inherent in tomographic binning, and we also show that the sensitivity 

O , to systematic errors in photometric redshift is much less. In conjunction with the fact 

■ ■ that the physics of lensing is very soundly based, the analysis here demonstrates that 

deep, wide-angle 3D weak lensing surveys are extremely promising for measuring dark 
energy properties. 
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1 INTRODUCTION 

Our knowledge of cosmology has advanced considerably in 
recent years. Led by detailed measurements of the microwave 
background radiation and large-scale structure, many of the 
most important cosmological parameters are now known 
with good accuracy. This advance has come about princi- 
pally through the all-sky maps of the microwave sky taken 
with the Wilkinson Microwave Anisotropy Probe (WMAP) 
(Bennett et al., 2003), supplemented by higher-resolution 
observations of the Arcminute Cosmology Bolometer Array 
Receiver (ACBAR) and Cosmic Background Imager (CBI) 
(Kuo et al., 2004; Pearson et al., 2003). When combined with 
large-scale structure information from the Anglo- Australian 

2 degree field galaxy redshift survey (2dFGRS) (CoUess et 
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al., 2001; Percival et al., 2001), the Lyman-a forest (Croft et 
al., 2002; Gnedin and Hamilton, 2002), and measurements of 
galaxy bias (Verde et al., 2002), the data establish the con- 
cordance model of an accelerating Universe dominated by 
dark energy and dark matter (Spergel et al., 2003; Spergel 
et al., 2006). The acceleration of the Universe is also appar- 
ent in observations of distant supernovae (e.g. Riess et al., 
2000). The determination of the density, baryon content and 
expansion rate of the Universe shifts the major unanswered 
questions in cosmology to the nature of the dark matter 
and dark energy. Dark energy in particular can be probed 
through its cosmological effects on the distance-redshift re- 
lation and the growth rate of structure. 

The question of the precise nature of the dark energy 
is a far-reaching one. We use the simplest phenomenological 
model of dark energy by parameterizing the equation of state 
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of the vacuum, 

w = p/{pc^), (1) 

where p is its energy-density and p is the dark- 
energy/vacuum-pressure. w = w{a) may vary with scale fac- 
tor. If it is found with high precision that w = — 1, then the 
dark energy cannot be distinguished with large-scale mea- 
surements from a modification to the gravity law along the 
lines suggested by Einstein with the cosmological constant. 
If, however, it can be established with a degree of certainty 
that w differs from —1 at any redshift, then it cannot be as- 
sociated with such a change to the gravity law, and is most 
naturally accounted for by a now field. This would be an ex- 
tremely important discovery, and the time-evolution of the 
field would be a useful constraint on models. Some possibil- 
ities exist in the literature, such as those proposed by Ratra 
and Peebles (1988), but none is a clear favourite candidate. 

Weak lensing is a very attractive proposition for study- 
ing dark energy, as it is sensitive to both of these effects, 
and, equally importantly, the physics of weak lensing is well 
understood. A key part of this is that it is sensitive to the 
distribution of matter in the Universe, regardless of its form. 
Furthermore, since weak lensing analysis can be done in a 
way which is either dependent on the distance-redshift rela- 
tion alone (see e.g. Taylor et al., 2006; Jain & Taylor, 2003) 
or on both the distance-redshift relation and the growth fac- 
tor (this paper), it can in principle distinguish between mod- 
ified gravity models and dark energy models. 

The main lesson of the field of microwave background 
astronomy is that with well-understood physics, robust re- 
sults can be obtained with high precision. Weak lensing ob- 
servations are, however, a technical challenge, as the imag- 
ing requirements are severe. Thus, it has only been in the 
last five years or so that the first measurements of cos- 
mic shear have appeared (Bacon, Rcfrcgier and Ellis, 2000; 
Kaiser, Wilson and Luppino, 2000; van Wacrbckc ct al., 
2000; Wittman et al., 2000). Weak lensing measurements 
to date have concentrated on obtaining the matter density 
parameter Urn and the amplitude of mass density fiuctua- 
tions (Hoekstra, Yee and Gladders, 2002; Jarvis et al., 2003; 
Rhodes et al., 2004; Hcymans et al., 2004; Hoekstra et al., 
2006; Semboloni et al., 2006). More ambitiously, weak lens- 
ing observations have started to put constraints on the equa- 
tion of state of dark energy (Jarvis et al., 2005; Semboloni at 
al., 2006). Theoretically, the prospects for determining dark 
energy properties (specifically its equation of state w) using 
weaJi lensing have been explored in a number of papers (e.g. 
Taylor et al., 2006; Hu and Tegmark, 1999; Huterer, 2002; 
Heavens, 2003; Rcfrcgier, 2003; Simon, King and Schneider, 
2004; Takada and Jain, 2004; Song and Knox, 2004; Ishak et 
al., 2004; Ishak, 2005). The prospects for determining w as a 
function of redshift z are markedly improved when 3D infor- 
mation on the individual lensed sources is available. Source 
distances could come from spectroscopic redshifts, but given 
the depth and the sky area required, they are more likely to 
be estimated from photometric redshifts. With 3D informa- 
tion, the lensing pattern can be analyzed in shells at different 
distances (e.g. Hu, 1999; Hu and Jain, 2004; Ishak, 2005), or 
by analyzing the shear pattern as a fully three-dimensional 
field (Heavens, 2003). It is the latter possibility which we 
investigate in this paper. The statistical properties of the 
shear pattern are influenced by many cosmological parame- 



ters, including w{z). In this paper we extend the analysis of 
Heavens (2003) to small-angle surveys as well as computing 
the expected marginal errors on w (and its evolution), using 
a Fisher matrix approach. We investigate issues of depth vs 
area, and the number of photometric bands which should be 
used, to determine the dark energy properties as accurately 
as possible. The main focus of the paper is in computing 
the expected statistical errors, but we do consider the im- 
pact of some systematics (Ishak et al., 2004; Bernstein, 2005; 
Huterer et al., 2005). 

The layout of the paper is as follows: in Section 2 we de- 
tail the transform method used and compute the covariance 
matrix of the transform coefficients; in Section 3 wc outline 
how the expected statistical errors on parameters are calcu- 
lated; in Section 4 we present the survey design and how we 
can scale to other surveys and in Section 5 we present an 
optimization of survey design and the parameter errors; in 
Section 6 we consider the synergy of 3D weak lensing with 
other dark energy probes and discuss future surveys and 
finally we give our conclusions in Section 7. 



2 METHOD 

2.1 Transformation of scalar and shear fields 

The observable quantities we use are the estimates of the 
shear field at locations in three dimensions. The estimates 
of the complex shear come from the shape and orientation 
of galaxies, where the radial distance is obtained approx- 
imately by using photometric redshift estimates obtained 
from observations through several or many filters. 

In a previous paper (Heavens, 2003) we introduced the 
idea of 3D weak lensing analysis in harmonic space as a 
statistical tool. In Castro, Heavens and Kitching (2005), we 
developed the subject formally and found the power spec- 
trum of 3D weak lensing shear. In this paper we consider the 
fiat-sky limit including the non-linear evolution of the power 
spectrum. We consider a transform of the 3D shear field in 
spin-weight spherical harmonics and spherical Bessel func- 
tions. This is a very natural expansion for the shear field, 
as the complex shear 7 is a spin-weight 2 object, as arc the 
spin-weight 2 spherical harmonics: under a local rotation of 
the coordinate system by angle ip, 7 changes to je^'^ . The 
spherical harmonic transform of a spin- weight s field s/(r) 
is defined here by 

Jem{k) = f rf'r»/(r) kjeikr) sY/^*{n) (2) 

where je{z) is a spherical Bcsscl function, sY™ a spin- weight 
s spherical harmonic, is a radial wavenumber, £ is a pos- 
itive integer, rn = —£,...£ and n represents the direction 
9, (p. For s = the spin-weight spherical harmonics are the 
usual spherical harmonics 17", and this is the appropriate 
spherical expansion of a scalar field. Note the presence here 
of a benign factor of k, to agree with the notation of Cas- 
tro, Heavens and Kitching (2005). The motivation for using 
spherical coordinates is manyfold: firstly the selection func- 
tion for a survey can often be separated into an angular 
(sky coverage) part and a radial component; secondly the 
errors in photometric redshifts introduce purely radial er- 
rors in the positions of the source galaxies; thirdly, in the 
Born approximation, the lensing effect is an integral effect 
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along the (radial) line of sight. The motivation in flat space 
for using products of spherical Bessel functions and spher- 
ical harmonics is that, as eigenfunctions of the Laplacian 
operator, it is easy to relate the expansion coefficients of the 
gravitational potential to those of the density field. Similar 
considerations led Heavens and Taylor (1995) (see also Fisher 
et al, 1994; Tadros et al., 1999; Percival et al., 2004) to ex- 
pand the large-scale structure of galaxies in spherical Bessel 
functions and spherical harmonics. Since cosmic shear de- 
pends on the gravitational potential, the use of this basis 
allows us to relate the expansion of the shear field to the 
expansion of the mass density field. The properties of the 
latter depend in a calculable way on cosmological parame- 
ters, so this opens up the possibility of using 3D weak shear 
to estimate these quantities. 

For surveys with large opening angles on the sky, a full 
expansion in spherical Bessel functions and spherical har- 
monics is the natural choice. Such an expansion is generally 
applicable, but for small-angle surveys whose signal is dom- 
inated by high-^ modes, the spherical harmonics are cum- 
bersome and their accurate computation can present prob- 
lems. For such surveys, we can approximate the spherical 
harmonics as sums of exponentials, as detailed in Appendix 
A of Santos ct al. (2003). 

For this paper we use the flat-sky expansion, which for 
a scalar (s = 0) field reads 



(3) 



where £ is a 2D angular wavenumber and k a radial 
wavenumber. In the spherical Bessel function, £ = \£\; £ is 
necessarily an integer, but we assume that £ ^ 1 so that 
enforcing integer ^ is a minor approximation. Note that we 
arc performing a full 3D expansion of the shear field and as- 
sume a flat Universe except where indicated. An alternative 
approach to include at least some 3D information is what is 
referred to as tomography, where the shear pattern of galax- 
ies is analyzed in shells, based on their photometric redshifts 
(Hu, 1999; Hu, 2002; Jain and Taylor, 2003; Takada and 
White, 2004). It is however evident that the binning process 
loses at least some information, and it is not necessary. 

The inverse transform in the flat-sky approximation is 



/(r) = 



TV J (2n] 



-dkkji{kr)e-xp{i£.e)f{k,£). (4) 



The coefficients of the expansion in the two systems are 
related by generalization of equation (A13) in Santos et al. 
(2003): 




(5) 



where the small survey is centred at the pole of the co- 
ordinate system, and the 2D transverse wavevector is £ = 
(£cos (j>e,, i sin (j)^). The covariances of the flat-sky coefficients 
are related to the power spectrum of / by 

(/(fc, £)r {k',£')) = (2^)^ Pfik) 5°(fe - k')5''{£ - £') (6) 

where is the Dirac delta function. 

Our plan is essentially to transform the components of 
the 3D shear field to produce a set of transform coefficients 



as a function of (k,£). These data will depend on cosmolog- 
ical parameters, and can be used in a likelihood analysis to 
constrain those parameters. 



2.1.1 Transformation of shear fields 

The weak lensing shear components we transform are 71 (r) 

and 72 (r), which are related to the lensing potential 0(r) 
through (e.g. Bartelmann and Schneider, 2001) 



71 = 2 ('^11 



4>22) ; 



72 



(7) 



where = d'^(f>/d9id9j. (f> itself is dependent on cosmologi- 
cal parameters through its relation to the mass density field 
(see section 2.4). We will return to this dependence later. 
For a large-area survey, it is a measure of the shears with 
respect to axes based on the spherical coordinate system, 
in which case the complex shear 7 s 71 -|- 172 is the second 
edth derivative of d>: 



7(r) 



;g g0(r) 



(8) 



(Castro, Heavens and Kitching, 2005). In the flat-sky limit, 
9 — > —{dx + idy), where the dx,y = d/d6x,y. Expanding 
the lensing potential in terms of spherical Bessel functions 
and exponential functions, as in equation (3), we see that 
it is natural to expand the complex shear field in terms of 
g 3 exp(-z£.0) = £^X^expi-i£.e), where 



X£^ 



(9) 



The £^ in the denominator is included for convenience, so the 
inverse transform kernel is just yj2/'Kk jt,{kr)X^ exp{i£.0g). 

2.1.2 Fiducial Cosmology 

An immediate issue to address is which radial coordinate to 
use in the spherical Bessel function. The observed quantities 
are the estimated redshifts of the sources, and we need to 
do two things: one is to translate these into radial distances; 
the second is to account for the error in the estimation of 
the redshifts. For the former, we choose a fiducial set of 
cosmological parameters, to define a transformation r^{zp) 
from the photometric redshift estimate Zp to a radial coordi- 
nate r". For this paper, we choose as the fiducial model the 
concordance model (Spergel et al., 2006) with fim = 0.27, 
= 0.04, = 0.73, 0-8 = 0.8, h = 0.71, wo = -1 and 
Wa = where the variables are the matter, baryon, vac- 
uum density parameters, Hubble constant in units of 100 
kms~'^Mpc~^ and the dark energy equation of state pa- 
rameters respectively. The equation of state of dark energy 
is modelled in terms of scale factor a by 



«;(a) = w -\- Wa(l — a) 



(10) 



(Chevallier and Polarski, 2001; Linder, 2003) where a{z) = 
(1 + z)^^ is the cosmic scale factor normalized to unity at 
the present epoch. 

We also include the scalar spectral index Us = 1 and 
its running a„ = 0. For the CMB Fisher calculations, see 
Section 5.1, we also include the tensor to scalar ratio r = 
0.01 and the optical depth to the surface of last scattering 
r = 0.09. 
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2.2 Transformation 

The lensing potential is defined everywhere, but we sample 
it only at the locations of galaxies, so it is natural to make 
a transformation of this point process, summing over galax- 
ies rather than integrating over space. Our estimate of the 
transform is thus defined as 



(11) 



where W{r) is au arbitrary weight function, and {rg,0g) are 
the coordinates of galaxy g. 

Note the appearance of two distances in the transform, 
r and (at each galaxy g): the main application of this 
study is to determine cosmological parameters, which affects 
the r{z) relation. The shear field is the shear field at the 
actual coordinate of the galaxy, and this depends on the 
true cosmological parameters, whereas the expansion (and 
weighting) is done with the fiducial model parameters. This 
distinction was neglected in Heavens (2003) and leads to 
an underestimate of the errors on the dark energy equation 
of state in that paper; the error estimates for the power 
spectrum in that paper are unaffected by this error. 

Writing the number density of source galaxies n(r) as 
the sum of a set of delta functions, we see that 

^{k,£) = ^ J d^rn{r)^{r)kje{kr'^)exp{-i£-0)W°, 

(12) 

where = W{r°). Note that in the high-^ hmit these are 
also the (minus) coefficients of the expansion of the conver- 
gence field K (Castro, Heavens and Kitching, 2005). This has 
an expectation value which is obtained by replacing rt(r) by 
the mean density of the source galaxies, n(r). Here we as- 
sume that selection effects are uniform across the survey so 
there is no angular dependence. Thus the 7 are estimators 
of 



7(fc,^) 



(frn{r)f{r)kje{kr°)exp{-ie.e)W°. (13) 



The estimates will differ because of the discrete nature of the 
galaxies, which leads to shot noise, the photometric redshift 
errors, and the source clustering. For deep surveys, and with 
a radial smoothing arising from the photometric redshifts, 
source clustering can be safely ignored. We include the ef- 
fects of photometric errors, but ignore uncertainties in the 
photometric redshift distribution. In terms of the observ- 
able photometric redshift distribution of sources (all-sky), 
n{r)d^r = nz{zp)dzp/4,%, we have 

^{k,tj = J dZpd^en,{zpMr)kMkr°) (14) 

exp{-i£ ■ e)W°. 



2.3 Photometric redshift errors 

Photometric redshifts lead to a smoothing of the distribution 
in the radial direction. If we denote by p{zp\z) the probabil- 
ity of the photometric redshift being Zp, given that the true 



redshift is z, the mean of the expansion coefficients will be 

l{k,£) = \j j dzdzpd^0p(zp\z)nz{zp) (15) 

7(r) k3e{kr°) exp{-it ■ e)W° . 

Note that p{zp\z) is arbitrary; it will generally have a dis- 
persion which depends on redshift, and can, if desired, in- 
clude broad wings to account for a small percentage of catas- 
trophic failures in the photometric redshift estimates. We 
assume a Gaussian, with a 2-dependent dispersion: 



p{zp\z) 



exp 



(Zp Z -\- 2^bias) 



2tTi(.z) 



(16) 



Zbias is a possible bias in the photometric redshift calibra- 
tion, the effect of this on dark energy parameters is discussed 
in Section 5.5.1. Strictly the shear is estimated at the actual 
radial coordinate of the galaxy, which may differ from r{z) 
because of peculiar velocities. We can safely ignore these, 
whose effect is small compared with current photometric 
redshift errors. 



2.4 Relationship of 7(fe,-£) to cosmological 
parameters 

The lensing potential is related to the peculiar gravi- 
tational potential $ by a radial line-of-sight integral (e.g. 
Bartelmann and Schneider, 2001): 

0(r) = 4/' dr'FK{ry)^{r'). (17) 
c Jo 

where 

FK{r,r') = {/K(r-r')/[/K(r)/K(r')]}, and /K(r)d^' is 
the dimensionless transverse comoving separation for points 
separated by an angle dtp. The Robertson- Walker metric 
may be written ds^ = c^dt^ - R^{t) [dr^ + f]i{r)dip'^], and 
fKij) takes the values sinr, r, sinhr for curvature values 
k = 1,0,-1. For a flat Universe FK(:r,r') = (1/r' - 1/r). 

The peculiar gravitational potential is related to the 
overdensity field S{r) = [p(r) — p] /p by the comoving Pois- 
son's equation 



= ^^'"-^0 5 



2a{t) 



(18) 



where Ctm is the present-day matter density parameter. Ho is 
the present Hubble constant and a{t) = R{t)/Ro = 1/(1-1-2;) 
is the scale factor. 

Note that 5 itself is not a homogeneous field, because it 
evolves with time, and hence with distance from the observer 
through the light travel time. We get around the subtleties 
of this by defining at each epoch a homogeneous field by re- 
ferring all field measurements to that time. Thus we can, for 
example, define a power spectrum which is time-dependent, 
and hence r-dependent. This may seem a little strange, since 
we have transformed from r space. The transforms of the ho- 
mogeneous fields will be denoted by 5{k,t;r) etc. 

For high I, the transforms of $ and 6 (referred to epoch 
t or equivalently r) are related simply by 



$(/c,^;r) 



\5(k,l;r) 



(19) 



2k'^a(t) 

Inserting these definitions in the equation for 7(A;,^),we find 
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the relationship between 7(fc, £) and the transform of S: 

OQ tt2 roG 

^{k,£) = -^^f- 1 dzdz,v{z^\z) (20) 

y d'^ekje{kr'')W° n,{zp)e^p{-ie.0) 

/ dr'a^^ (r) FK(r,r) 
Jo 

Integration over gives (27r)^(5^(£ — £'), so 

7(fe,^) = — / dzdzpp{zp\z)n^{zp) (21) 

47r^c^ Jo 

je{kr°)W° f dr'FK{r,r'){l + z') 
Jo 

J dk' k'je{k'r')S{k',£;r'). 

This is a fundamental result of this paper. It establishes 

the connection between the (observable) 3D shear transform 
coefficients, and the underlying matter density fluctuations, 
whose properties are calculable from theory. 

2.5 Covariance matrix of 'y{k,£) 

The signal part of the covariance matrix of the 'Y{k,£) is 
obtained from equation (21). For the covariance of the over- 
density field coefficients, it is algebraically convenient to use 
the geometric mean of the power spectra Ps, rather than the 
power spectrum evaluated at epochs corresponding to r or 
r' . Both of these could also be justified; note also that Ps 
does not depend on £ (Castro, Heavens and Kitching, 2005) . 

{5ik,£;r)6'{k',£';r')) ~ i2nf Psik;r)Ps(k';r') (22) 

5'^ {k-k')6'^ {£-£'). 

The covariance matrix for the shear expansion coefficients 
is then 

(7(fc, £)7*(fc', i'))s = Qi {k, k') 5°{£ - £') (23) 
where Qi {k, k') can be written as 

Qdk^ k ) = J^2°4 ^ J ^ Ge{k, k)Ge{k', k) (24) 
where 

Ge{k,k) = J dzdZpnz{zp)W{zp)p{zp\z)Ue{r,k)je{kr°) 

(25) 

and 

Ue{r,k)= [ df ^''^X^ VPi{k;f)je{kr). (26) 
Jo 

where r = r{z) etc. Equation (23) is the second important 
result of the paper. 

2.6 Shot noise 

The shot noise can be calculated by making the usual as- 
sumption that the galaxies are a Poisson sampling of an 
underlying smooth field (see e.g. Peebles, 1980). In practice 
we consider estimators of the transforms of the individual 



components of the shear, 7q,; a = 1,2. In the normal way 
for a point process, these may be written as sums over small 
cells c, each of which contains ric = or 1 galaxy: 

^a{k,£) = \ - nc'yac{rc)kje{krc)exp{-i£-0c)W°. 

V ^ cellTc 

(27) 

The variance of this involves a double sum over cells, and 
the averaging over cells c and d, {ricnd) contains shot noise 
terms when c = d, in which case (ric) = («c), and the shot 
noise reduces to a single sum, or an integral when we move 
back to a continuum description. Using the fact that the 
variance of the shear estimate for a single galaxy is com- 
pletely dominated by the variance in the intrinsic ellipticity 
of the galaxy, cTg, rather than by lensing, 

M) = (28) 

where 5^ is a Kronecker delta function, and Ue ~ 0.3 (Brown 
et al., 2003), we find an expression for the shot noise as 

{%{k,£}%{k',£'})sN = <tI j dzn,{z)kje{kr°) (29) 

k'Mk'r°)W\z)dSg6''i£-£'). 



3 ESTIMATION OF COSMOLOGICAL 
PARAMETERS 

Cosmological parameters infiuence the shear transforms in 
a number of ways: the matter power spectrum Ps{k;t) is 
dependent on Qrn, h and the linear amplitude erg. The linear 
power spectrum is dependent on the growth rate, which also 
has some sensitivity to the vacuum energy equation of state 
parameter w{z). w also affects the r{z) relation and hence 
the angular diameter distance fK[r{z)]. These parameters 
({Oa}) may be estimated from the data using likelihood 
methods. Assuming uniform priors for the parameters, the 
maximum a posteriori probability for the parameters is given 
by the maximum likelihood solution. We use a Gaussian 
likelihood 

2\nL{g\{9c}) = constant - det(C) - g • C"^ • g (30) 

where C = S + N is the covariance matrix, given by sig- 
nal and noise terms equations (23) and (29). Note that the 
average values of 7(fc, £) is zero, so the information on the 
parameters comes from the dependence of the signal part of 
the covariance matrix C. i.e. we adjust the parameters until 
the covariance of the model matches that of the data. This 
was the approach of Heavens and Taylor (1995); Ballinger, 
Heavens and Taylor (1995); Tadros et al. (1999); Percival et 
al. (2004) in analysis of large-scale galaxy data. For many 
surveys, many useful modes of the shear transform have con- 
tributions from wavenumbers where the power spectrum is 
quite nonlinear. The use of a Gaussian likelihood thus needs 
to be justified by comparison with simulated data; this will 
be explored in a later paper, and it is possible that a differ- 
ent likelihood function may be necessary in the non-linear 
regime. 
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3.1 Expected errors on cosmological parameters - 
the Fisher matrix 

The expected errors on the parameters can bo estimated 
with the Fisher information matrix (Jungman et al., 1996; 
Tegmark, Taylor and Heavens, 1997). This has the great 

advantage that different observational strategics can bo an- 
alyzed and this can be very valuable for experimental de- 
sign. The Fisher matrix gives the best errors to expect, and 
should be accurate if the hkelihood surface near the peak is 
adequately approximated by a multivariate Gaussian. 

The Fisher matrix is the expectation value of the second 
derivative of the InL with respect to the parameters Oa'- 



(31) 



and the marginal error on parameter Qq, is ^ {F~^")aa- If 
the means of the data are fixed, the Fisher matrix can be 
calculated from the covariance matrix and its derivatives 
(Tegmark, Taylor and Heavens, 1997) by 



= lli-ace [C-'C,^C-'C,g] . 



(32) 



For a square patch of sky, the Fourier transform leads to 
uncorrelated modes, provided the modes are separated by 
27r/L where L is the side of the square in radians, and the 
Fisher matrix is simply the sum of the Fisher matrices of 
each £ mode: 



Fc0 = IJ2 Trace [icT'c%icY'c: 



131 



(33) 



where is the covariance matrix for a given £ mode. Wc 
compute numerically from the signal and noise parts 
equations (23) and (29), for given nz{z), photometric red- 
shift error distribution, cosmology and survey area, which 
governs the separation of uncorrelated i modes. 



4 SURVEY DESIGN FORMALISM 

In this Section we discuss the survey design factors. We start 
by detailing the assumptions of the survey design, and dis- 
cuss some details of the Fisher matrix calculation. Possible 
future weak Icnsing surveys and their effectiveness are dis- 
cussed in the Section 5. 



4.1 Survey parameters 

In assigning survey parameters wc follow the formalism de- 
tailed in Taylor et al. (2006). We assume that the redshift 
distribution for a typical magnitude-limited survey is of the 
form 



n{z) oc z exp 



(34) 



where zq = 2^/1.412, and Zm is the median redshift of the 
survey (e.g. Baugh and Efstathiou, 1993). The number den- 
sity of useable sources with photometric redshift and shape 
estimates is taken to scale as 



no = 30z^ per square arcminute. 



(35) 



This was estimated from the COMBO- 17 survey. We take a 
maximum redshift of Zmax = 1-5 for ground based surveys. 




Redshift 

Figure 1. The photometric errors as a function of redshift z for 
a 5-band (upper sohd) and a 17-band (lower solid) optical survey. 
The dashed line has (Tz{z) = 0.05(1 -|- z) for compaxison. 



This is due to the difficulty of measuring galaxy shape, be- 
cause of the decrease in a galax;y's apparent size with in- 
creasing redshift, coupled with the seeing limit. 

We assume that the photometric redshift errors are 
Gaussian, with a dispersion given by cr{z) = <to(1 + 
z) f(m,, z). where m. is the apparent magnitude of the galaxy, 
the function /(m, z) is given in Taylor et al. (2006). We in- 
tegrate over a Schechter function to get the average error 
as a function of z. The error distribution is shown for a 5- 
band optical survey and a 17-band optical survey in Figure 
1. The assumption of Gaussianity of p(0p|z) can easily be re- 
laxed: outliers can, for example, be included, we investigate 
this in Section 5.5.2. We have extrapolated these formula 
to 2: = 1.5, though this extrapolation may be optimistic as 
photometric redshift estimates can increases dramatically at 
z ~ 1 if IR data is not available. 

The variables which can be varied are the area A and 
depth of the survey (zm), and the number of bands. These 
scale with the number of nights observing T . the telescope 
diameter D and the field-of-view F as (see Taylor et al., 
2006) 



T oLztAD-^F- 



(36) 



We consider as our default survey a 4 metre telescope with 

a 2 square degree field-of-view which could observe an area 
of 10,000 square degrees to Zra ~ 0.7 with 5-bands in 600 
nights of observing, this could be achievable with surveys 
such as darkCAM (Taylor, 2005; conference proceedings of 
Probing the Dark Universe with Subaru and Gemini) or the 
Dark Energy Survey (Wester, 2005). 

We compute the nonlinear power spectrum using the fit- 
ting formulae of Smith et al. (2003), based on linear growth 
rates given by Linder and Jenkins (2003). In order to avoid 
the high-fe regime where the formulae may be unreliable, 
or where baryonic effects might alter the power spectrum 
(fe > lOhMpc'^; White, 2004; Zhan & Knox, 2004), we do 
not analyse modes with k > 1.5 Mpc~^ . Note that the non- 
local nature of lensing does mix modes to some degree, but 
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I mode 

Figure 2. The contributions to tlic Fislicr matrix element Fyjyj 
from different ranges of £ for a 10,000 square degree survey with 
Zm = 0-70 and 5 bands, w = wo in this Figure for clarity. 



Lensing Survey 

Area/sq degrees Zm z^ax ^'^Bands 

10,000 0.70 1.5 5 



Planck 

Band/GHz 6>beam (Tt/IQ-^ ctp/IO"' 

44 23' 2.4 3.4 

70 14' 3.6 5.1 

143 8.0' 2.0 3.7 

217 5.5' 4.3 8.9 



WFMOS 

Area/sq degrees zy^in femax /hMpc~^ Bias 

2000 1.0 0.15 1.25 

300 3.0 0.15 1.25 



these modes are sufficiently far from the uncertain highly- 
nonhnear regime that this is not a concern (Castro, Heavens 
and Kitching, 2005). In any case the results are very insen- 
sitive to the radial k limit, since the photometric redshift 
errors suprcss radial power at much lower k. We include an- 
gular modes as small as each survey will allow, and analyse 
up to ^max = 5000. We note that the intracluster medium 
might affect the power spectrum on the level of a few percent 
for 1000 <£< 3000 (Zhan & Knox, 2004). These modes will 
still contain useful information, but a more detailed analy- 
sis might be ncssacary when the method is applied. To help 
asses the extent of any modification to the expected accu- 
racy, we also quote results for a more conservative limit of 
^max = 2000, this increases the predicted marginalised er- 
rors by approximately 0.01. The flat-sky approximation will 
break down for the low-^ modes, but there is little power 
there in any case (Figure 2). 

We allow for a Universe with the following parameters: 
0,m, ^iv, h, ug, Qb, Too, Wa, Us and a„. as represents the 
amplitude of the perturbations, Us the scalar spectral index 
and its running q„, and we parameterize w{a) by equation 
(10). Note that our assumption of this form is not critical; 
theoretical models with arbitrary w{a) can be analyzed. We 
clioosc> this form to investigate the sensitivity of our results 
on Wo and Wa to time-dependence. 



5 PARAMETER FORECASTS AND 

OPTIMIZATION FOR A WIDE-FIELD 
LENSING SURVEY 

Having introduced the method and the survey design formal- 
ism this Section will investigating optimizing a weak lensing 
survey so that the marginal errors on the dark energy pa- 
rameters can be minimized. We will explore the variation in 
the marginal error on wo with changes in the median depth, 
varying the area to preserve the total observation time, and 
the redshift error. 

© 0000 RAS, MNRAS 000, 000-000 



SNAP 


^max 




NsNIa 


Cm 


1.5 


17 


2000 


0.15 



Table 1. The main default values parameterising the Lensing, 
CMB, BAO and SNIa experiments considered in this paper. 



5.1 Combining with other dark energy 
experiments 

There are a number of alternative ways to place constraints 
on the dark energy equation of state parameters. Prominent 
among these are the CMB, Baryon Acoustic Oscillations 
(BAO) in the galaxy power spectrum and the supernova 
Type la (SNIa) Hubble diagram. Each of these methods can 
place constraints on the dark energy equation of state al- 
though they all suffer from degeneracies between wo and 
Wa- In Section 6 we show how by combining the different 
methods the parameter degeneracies can be lifted allowing 
for improved accuracies on wo and Wa- 

Using the Fisher matrix methods outlined in detail in 
Taylor ct al. (2006) we calculate predicted Fisher matri- 
ces and parameter constraints for the following surveys. In 
all the Fisher matrix calculations we use an 11 parameter 
cosmological set {f^m, i^v, h, ag, i^b, wo, Wa, ris, r, a„, 
r = T/S), with default values (0.27, 0.73, 0.71, 0.8, 0.04, 
-1.0,0.0, 1.0, 0.09, 0.0, 0.01). For a summary of the main 
assumptions that went into each of the Fisher matrix calcu- 
lations see Table 1. 

We consider a 4-year WMAP experiment and a 14- 
month Planck experiment (Lamarre et al., 2003) which will 
be contemporary with the type of wide field photometric 
surveys being considered for 3D weak lensing. We com- 
puted the CMB covariance matrices by using CMBFAST 
(Seljak and Zaldarriaga, 1996), the survey parameters were 
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median 

Figure 3. The variation in the marginal error on wq as the me- 
dian redshift of the survey varies for a 600 night survey on a 4 
metre elass telescope, including a 14-month Planck prior. Note we 
assume shapes are note measurable beyond .Smax — 1.5. The solid 
line is for a 5 band survey, the dashed line for a 9 band survey 
and the dot-dashed line for a 17 band survey. 



taken from Hu (2002). Wc include polarisation but do not 
marginalize over the calibration error. Also we do not con- 
sider the Integrated Sachs- Wolfe (ISW) effect directly via 
cross-correlating with galaxy surveys, although this will also 
provide an interesting dark energy constraint. 

For a BAO experiment we consider the Wide Field 
Multi Object Spectrograph (WFMOS; Bassett et al., 2005) 
following the methodology outlined in Seo & Eisenstein 
(2002), Blake & Glazebrook (2003) and Wang (2006). 

We use the method outlined in Ishak (2005) and Yeche 
ct al. (2006) to calculate a SNAP (Aldering, 2005) SNIa 
Fisher matrix. The effective magnitude uncertainty takes 
into account luminosity evolution, gravitational lensing, 
dust and the eflfect of peculiar velocities. 



5.2 Combining with the CMB alone 

To help to lift degeneracies between cosmological parame- 
ters, and to retain realism in our predictions for a wide field 
photometric survey we consider a 14-month Planck CMB 
experiment. This CMB experiment will place constraints on 
Wo and Wa, mainly through the large scale ISW effect, al- 
though there is a strong degeneracy between the dark energy 
parameters. Combining with 3D weak lensing helps to lift 
this degeneracy. 

For the default survey, we show in Figure 4 the Fisher 
matrix elements marginalized over all other parameters. The 
darkest areas show a 14-month Planck prior. The light gray 
ellipses show the two-parameter, 1-a errors for the parame- 
ters plotted, and the white (central) ellipses show the combi- 
nation. The marginal errors on wo and Wa are Awo = 0.108 
and Awa = 0.395 respectively, a factor of 5 improvement 
over the 14-month Planck constraints alone which could con- 
strain Wo and Wa to Awo = 0.502 and Awa = 1.86. 

3D weak lensing improves constraints on all the CMB 



P n TP> m fit f^r 


P Ian r'K i~irt 1 


T.finGTncr r»nl\7 


l! r» TYi r» in <3 H 
V—"-! Ill 1 lie u 


n 


0058 


0500 


0025 


n 


0.0024 


0.0795 


0.0015 


h 


008S 


0.0321 


0.0051 




0.1002 


0.0705 


0.0073 


Oi, 


0.0011 


0.3707 


0.0007 


Wq 


0.5015 


0.2843 


0.1086 


Wa 


1.8618 


1.1792 


0.3947 


Us 


0.0034 


0.3852 


0.0031 


O-n 


0.0062 


0.0576 


0.0045 


T 


0.0079 




0.0077 


r 


0.0208 




0.0203 



Table 2. Improvements on CMB Planck one parameter l-u, con- 
straints by adding 3D weak lensing from a 10,000 square degree 
lensing survey to a median depth of Zm = 0.7. 



cosmological parameters, in particular as whose constraint 
is improved by a factor of 14. It is already well known that 
weak lensing can tightly constrain the {as, fim) plane, using 
standard cosmic shear techniques (see Brown et al., 2003; 
Semboloni et al., 2006), 3D weak lensing constrains as in 
the same way by measuring the overall normalisation of the 
matter power spectrum. 3D weak lensing provides indirect 
(and slight) improvements on the constraints for the ten- 
sor to scalar ratio r and the optical depth to last scatter- 
ing T through the intersection of the 3D weak lensing multi 
parameter ellipsoid with the CMB's multi parameter con- 
straint. Table 2 shows the 14- month Planck constraints and 
the new combined constraints once 3D weak lensing is in- 
cluded. Note that results arc presented for Universes which 
are not necessarily flat. In non-Hat geometries the spherical 
Bessel functions je{kr) should be replaced by ultra-spherical 
Bessel functions $^(j/). For the case considered here £ ^ 1 
and k ^ (curvature scale) then ^yiy) je{kr) (Abbott 
and Schaefer, 1986; Zaladarriaga and Seljak, 2000). The ex- 
pansion used is not ideal for non-flat Universes but should 
be an adequate approximation given current constraints on 
flatness. 



5.3 Survey Optimization 

For a given observing time, there will be an optimum depth 
of survey to minimize the statistical error on wo (or Wa)- A 
very wide, shallow survey will yield poor cosmological con- 
straints since the lensing signal is so small, whereas a very 
deep survey will also yield poor cosmological constraints be- 
cause very little area can be covered, and the cosmic vari- 
ance will be large. In addition, the distant galaxies will have 
shapes which are difficult to measure at high redshift. Here 
we explore the optimum median redshift using equation (36) 
keeping the time fixed, so that the area of the survey scales 
with median redshift as A cx . The results are for 600 
nights on a 4m survey telescope with a 2 square degree field- 
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Figure 4. Expected marginal errors on cosmological parameters from Planck (dark), 3D weaJc lensing survey (light gray) and the 
combination (central, white). The survey covers 10, 000 square degrees to a median depth of 0.7 in 5 bands. Ellipses show the 1-a errors 
for two parameters (68.3% confidence regions), marginalized over all other parameters. 



of-view, where Zm = 0.70 corresponds to A = 10, 000 square 
degrees. The results arc shown in Figure 3. 

The optimal median redshift for a 5 band optical sur- 
vey is Zm = 1.0 with an area oi A = 2400 square degrees. 
The error is poor below z,n ^ 1.0 because the lensing signal 
becomes weaker, as the number of Icnscd background galax- 
ies at high redshift decreases. For median redshifts above 
«m ~ 1.0 the error is also poor as shot noise begins to dom- 
inate and the areal coverage becomes too small. We also 
investigate the effect of using 9 and 17 optical bands, keep- 
ing the median redshift and area fixed. Wc find that the 
marginal error does not substantially decrease. This is due 
to the combination of the lensing and CMB, the intersec- 
tion of the ellipses in parameter space, remaining similar 
even though the lensing marginal errors on their own de- 
crease with increasing bands. The extra bands might well 
be useful in the identification of outliers which we discuss in 
Section 5.5.2. 

The optimal median redshift increases with the number 
of optical bands, a 17 band optical survey has an optimal 
median redshift of = 1-1 with an area of j4 = 1640 square 



degrees. Higher-ft modes are accessible due to the reduced 
damping effect of the photometric redshift error. In order 
to utilise these modes, the shot noise needs to be reduced, 
so the optimisation favours a slightly deeper survey with 
higher number density. This effect increases as the redshift 
uncertainty decreases. 



Taylor et al. (2006) find that the optimal survey design 
for the geometric ratio method is Zm = 0.70 and A = 10, 000 

square degrees for a 5 band photometric redshift survey on 
a 4 metre telescope with a 2 square degree field-of-view. So 
that the results presented here can be directly comparable 
with their results we shall adopt a fiducial survey design 
of Zm = 0.70 and A = 10, 000 square degrees in 5 bands 
from hereon; although it should be noted that given 600 
nights on such an instrument that an optimal survey design 
could improve the marginal errors on wo and Wa using this 
method. 
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Figure 5. The variation in the lensing-only marginal error on wq, 
as the weighting scheme is varied. The solid line is for W{z) = z", 
the dashed line is for W{z) = (1 + z)"; a is varied. 



Figure 6. Marginal error on wq for different photometric redshift 
errors parameterized by 1x^(2) = cro(l + 2). These results include 
a 14-month Planck prior. 



5.4 Weighting the data 

We might expect that the cosmological parameter con- 
straints could be improved if the more distant galaxies are 
given higher weight. We have not attempted a formal opti- 
mization, but show some results of experimental weighting 
schemes. We consider two weighting schemes W{z) = z" 
and W{z) = (1 -|- z)". Figure 5 shows how the error on 
Wo changes as the weighting scheme is changed. The figure 
shows the Icnsing-only marginal combined with a 14-month 
Planck prior on all parameters. We show that the marginal 
error on wo is fairly insensitive to the weighting scheme em- 
ployed. Furthermore using equal weighting is in fact the op- 
timal strategy for a weighting functional form of this kind. 
This shows that the increase in the shot noise through the 
weighting of high redshift galaxies counteracts any improve- 
ment in the lensing signal, used to constrain wq. 




5.5 Optical and Infrared surveys 

By combining a 5 band optical survey with, for example, 
a 4 band infrared survey the photometric redshift accuracy 
can decrease. Strategies such as this have the potential to 
be employed on future wide field surveys in an effort to im- 
prove cosmological parameter constraint. Here we parame- 
terize the redshift error using 



(Tziz) = (To(l +«) 



(37) 



where a 5 band optical survey can be approximately repre- 
sented, see Figure 1, by ao — 0.05. (Jq — 0.01 corresponds 
to a 9 band survey comprising of 5 optical and 4 infrared 
bands (Wolf, private communication). Note the distinction 
between this and the 9 band optical (no infrared) survey 
considered in Section 5.3. 

Figure 6 shows how the marginal error on wo varies 
with ao, here we fix the survey design to be Zm = 0.7 and 
A = 10, 000 square degrees. We find that the marginal er- 



Figure 7. Marginal error on wq for different photometric redshift 
errors parameterized by <Tz{z) = ito(1 -|- z) for lensing (solid line) 
and a BAO experiment (dashed line) from survey of Zm = 0.70 
and A = 10, 000 square degrees. These results include a 14^month 
Planck prior. 



ror on wo varies slowly between 0.01 < ao < 0.1 and im- 
proves rapidly for ao < 0.01. This turn-over corresponds 
to the point where the lensing pivot redshift error, see Sec- 
tion 5.7.1, becomes comparable to the CMB pivot redshift 
error. So, the marginal error of the combined constraint im- 
proves at a faster rate (as ctq decreases) after this point 
since the lensing constraint is improving wo and Wa and 
lifting tiie CMB degeneracies further. Figure 7 shows the 
marginal error from both lensing and BAO using the photo- 
metric redshifts from our fiducial survey. The treatment of 
a BAO experiment in a photometric redshift survey is dis- 
cussed in Taylor et al. (2006) . We find that BAO constraints 
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do not significantly improve the 14-month Planck prior for 
uo > 0.02, and that for a 5 band survey lensing provides 
much tighter constraints on the dark energy parameters. The 
dark energy constraints from lensing are less effected by the 
photometric damping of the radial wavcrmmbcr due to poor 
photometric errors than the BAO since lensing also gains 
dark energy information from geometric factor via the lens- 
ing effect. The BAO methodology on the other hand relies 
on a good measurement of the power spectrum which is re- 
stricted to low wavenumbers due to the photometric rcdshift 
errors. At low ao the two methods provide complimentary 
constraints on wq. 



only a single bias parameter, rather than the more complex 
behaviour allowed in Ma et al. (2005). However, it is the 
same model as used in our companion paper (Taylor et al., 
2006) which shows more sensitivity to Zbias. 

Following the procedure outlined in Taylor ct al. 
(2006) we also investigated the effect of an offset in 
the variance of the photometric redshift errors (Tz{z) — » 
y/ {criiz) + Aa'i{z)). We find that this effect is negligible for 
this analysis, so that the total bias due to photometric red- 
shift errors is only dependent on the bias in the offset of 
the mean. We will explore fully marginalizing over nuisance 
parameters in a full Fisher analysis elsewhere. 



5. 5. 1 Bias in the photometric redshifts 

As well as investigating the effect of varying the absolute 
values of the photometric errors, we would also like to in- 
vestigate how a bias in the photometric redshift calibration 
effects the dark energy parameter estimation. Knox, Scocci- 
marro and Dodelson (1998), Kim et al. (2004) and Taylor et 
al. (2006) show how the bias S-ipj on a fixed model param- 
eter tpj is related to the marginal error on a (cosmological) 
parameter 69i by 



(38) 



F^* is the cosmological parameter Fisher matrix and F^'^ is 
a pseudo-Fisher matrix between the measured and assumed 
parameters, when considering one parameter this is a col- 
umn matrix. 

We assume that there is some bias 2:bias in the mean 
of the photometric redshifts in a given survey (see equation 
(16)) due to poor calibration of photometric redshifts with 
a spectroscopic training set. Marginalizing over all other pa- 
rameters we find that 



5wo = —ChiaaSzbia 



(39) 



where Cbias is some constant. Following the arguments in 
Taylor et al. (2006) the number of galaxies requiring spec- 
troscopic redshifts is 



CbiasC(«) 



6wo 



(40) 



where the bias on wo is half the error 6wo = O.SAwq. 
We have found that Cbias ~ 1.2 for 3D weak lensing. If 
a{z) « 0.1 and we require Awo « 0.01 the number of spec- 



troscopic redshift required is Ng^ 



6 X 10^. This number 



is easily achievable using the current generation of spec- 
trometers. The number of required spectroscopic redshifts 
is significantly smaller than the large number required for 
the geometric ratio test, for which Cbias ~ 9.0 (e.g. Taylor 
et al., 2006) and tomographic methods (e.g. Hu and Jain, 
2004) , we attribute this difference to the binning procedure 
required in these methods. In binning the data any offset in 
the redshift estimation of a galaxy will create a discrepancy 
between the estimated and actual number of galaxies in a 
bin, and any derived quantities gained from them, for ex- 
ample the tangential shear behind a cluster. In this analysis 
any systematic offset in the galaxy population does not af- 
fect any derived quantities in this way but rather the whole 
shear field is offset in redshift, and galaxies are simply given 
a slightly increased weighting via je. Note that we consider 



5.5.2 The effect of outliers 

In every photometric redshift survey there will be some 

galaxies in the sample for which an accurate photometric 
redshift cannot be assigned. To investigate the effect of these 
outliers on the dark energy parameter estimation we con- 
sider two galaxy populations, one with the original photo- 
metric redshift errors, see Section 4.1, and a second popu- 
lation with a^^{z) — 0.5. We show the results in Figure 8, 
A'''^ is the proportion of the total galaxy population with 
crf^ = 0.5. We have considered three ways in which the out- 
lying population could be dealt with, and how the effect of 
each of these methods varies with the proportion of the total 
population of outlying galaxies. 

A population of outliers can either be discarded from 
the analysis completely or used in some way. The dashed 
line in Figure 8 shows the effect of discarding the sample, 
so that the surface number density of galaxies is decreased 
by no — > no(l — A''^), but the photometric redshift error 
remains the same. To use the outliers either they can be 
treated as a separate population (solid line in Figure 8) or 
can be incorporated into a single population (dot-dashed line 
in Figure 8) in which case the overall photometric redshift 

is degraded to ct,{z) ^J {1 - AP^)[a4z)]'' + AP^[af{z)]^ 
(see Blake & Bridle, 2005) where (yz{z) is the original pho- 
tometric redshift error. The effect of having outliers in the 
sample increases the marginal error on wo regardless of how 
they axe treated, though the method is relatively insensi- 
tive to this effect. As expected using the outlying galaxies, 
and treating them as a separate population, increases the 
marginal error less than discarding the galaxies completely. 
By incorporating the galaxies into a single population the 
redshift error is degraded to such a degree that for a low 
proportion of outliers it is optimal to discard them, note the 
signal-to-noise for 3D weak lensing is proportional to no. For 
a high proportion it is optimal to include them somehow ei- 
ther into a single population or as two separate populations. 



5.6 Scaling to other Surveys 

Using the results presented in Sections 5.3 and 5.5 we can 
provide scaling relations in a similar way to Taylor et al. 
(2006). To scale these results to other weak lensing surveys, 
equation (36) should be used with a time calibration i.e. 



To \zmoJ \AoJ \DoJ \fovoj 



(41) 
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Figure 8. The effect of outliers with (7^(z) = 0.5 on the marginal 
error in uJo as a function of the proportion of outliers in the survey 
AP'^. The survey is a 5 band survey of 10,000 square degrees to a 
median redshift of Zm = 0.7, with a 14-month Planck prior. The 

solid line shows the effect of treating the outliers as a separate 
population; the dashed line shows the effect of discarding the 
outliers; the dot-dashed line shows the effect of incorporating the 
outliers into a single gala:xy population. 





CFHT 


COMBO-17 


D(m) 


3.6 


2.2 


fov (sq deg.) 


3 


1 


N (bands) 


5 


17 


Zm 


1.17 


0.7 


Area (sq.deg.) 


170 


1 


T (nights) 


500 


6 



Table 3. Default survey parameters for the 5-band CFHT Legacy 
Survey and the 17-band COMBO-17 survey. 



The subscript refers to the parameters time, median red- 
shift and area of a survey on a telescope with certain diam- 
eter and field-of-view. The scaling applies between surveys 
with equal number of bands; for 5 bands the Canada- France- 
Hawaii Telescope Legacy Survey (CFHTLS) can be used, 
while for 17 bands COMBO-17 can be used. For other num- 
bers of bands it can bo naively assumed that the time for a 
given survey scales proportionally with the number of bands 
so that To ToNto/Nf, where Nf, is the number of bands in 
the survey. For a flexible survey design the optimal median 
redshift of «m ~ 1-0 is approximately insensitive to the num- 
ber of bands, when combined with a Planck prior (sec Figure 
3). If the number of bands is 5, 9 or 17 the appropriate line 
in Figure 3 then scales proportionally up (and down) with 
decreased (or increased) areal coverage from 2400 square de- 
grees, for a 5 band survey i.e. Awo(A) = (0.093) (A/2400) 
If the number of bands is not shown in Figure 3 then Fig- 
ure 6 can be used to find the minimum of the appropriate 
vs. 2^771 line (at Zm — 0.7). This can then be scaled for 
a differing areal coverage as before. 
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z 

Figure 9. Marginal error on w{z) combined with a 14-month 
Planck prior. The highest accuracy is achieved at the pivot red- 
shift of z = 0.373 with an error of Aw{z = 0.373) = 0.0175. 

For a fixed survey of area A and median redshift Zm the 
error in Figure 3 for a given median redshift Awo(^m) can be 
calculated using Awo{A) = Awo{zm){A/2400)-^. In scaling 
between bands a similar interpolation between Figure 3 and 
Figure 6 can be performed. 

5.7 Constraining w(z) at higher redshifts 
5.7.1 Pivot redshifts 

The parametrization used for the dark energy equation of 
state encodes information on both the present day value 
of w and its redshift evolution. By placing constraints on 
both Wo and Wa a region in the w{z) vs. redshift coordinate 
system is constrained. Through the anti-correlation of wo 
and Wa this constraint is minimized at a 'pivot redshift', the 
minimal error at this redshift being the pivot redshift error. 
The pivot redshift is schematically defined by the angle of 
the ellipse in the (wo,Wa) plane, the error at this redshift 
being the width of the semi-minor axis of the ellipse. 

Figure 9 shows the constraint on w{z) = wo + 'Wa{z/{l + 
z)) as a function of redshift for our fiducial weak lensing 
survey. The highest accuracy on w{z) occurs at the pivot 
redshift of ^ = 0.373 with Aw{z = 0.373) = 0.0175. The 
pivot redshift of the CMB alone is 2 = 0.368 with an error 
of Aw{z = 0.368) — 0.0350, the pivot redshift in this case 
is determined by the redshift at which the dark energy den- 
sity begins to dominate over the matter density. The pivot 
redshift of 3D weak lensing alone is z = 0.208 with an error 
of Aw{z = 0.208) = 0.2018 this is determined both by the 
redshift at which dark energy becomes dominant and the 
redshift at which the lensing signal is maximized. 

Figure 10 shows how the error on w{z) varies both with 
redshift and the median redshift of the survey. The line in 
Figure 9 can be found by tracing the Zm = 0.70 line in Fig- 
ure 10, the 5 band line in Figure 3 can be found by tracing 
along the x-axis {w{z) = w(0) = wo) in Figure 10. There is 
little sensitivity to the pivot redshift or the pivot redshift er- 
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median 

Figure 10. Marginal error on w{z) combined with a 14-month 
Planck prior as a function of median redshift. The contours arc 
lines of equal marginal error the values of the contour given on 
the line. 



ror on the survey design, this is due to the intersection of the 
3D weak lensing constraint with the 14-month Planck con- 
straint remaining the same. This occurs because the pivot 
redshift of 3D weak lensing is a property of the cosmological 
dependence of the method not the survey design parameters, 
so that despite the marginal errors on wo and Wa varying 
with the median redshift the orientation of the lensing el- 
lipse, and hence its intersection with the 14-month Planck 
ellipse, remains the same. 



5.7.2 Figure of Merit 

A 'figure of merit' has recently been introduced by the Dark 
Energy Task Force (DETF) (2006) (also see Linder, 2006; 
Taylor et al., 2006) which represents the area of the decor- 
related ellipse constrained by a survey at the pivot redshift 
and can be written 

Aw(^;pivot) * Awa. (42) 

In reference to Figure 9 the figure of merit quantifies both 
the minimal pivot redshift error and the redshift range over 
which the error in w{z) is small; a wide and deep curve in 
Figure 9 would have a small figure of merit. Figure 11 shows 
how the figure of merit for a survey consisting of 600 nights 
on a 4 metro telescope with a 2 square degree ficld-of-view 
varies with the median redshift of the survey. It can be seen 
that for a 5, 9 or 17 band optical survey the optimal median 
redshifts are the same as when optimizing for wo alone, see 
Figure 3. It can be seen that the optimal median redshift 
in Figure 11 coincides with the widest point of the inner 
contour in Figure 10. The Figure of merit is shown for all 
considered experiments in Table 4. 




median 

Figure 11. The figure of merit (product of errors in ^(zpivot) ^■iid 
Wa) as a function of median redshift for a Zm = 0.70, A = 10, 000 
square degree survey including a 14-month Planck prior. Errors in 

■wo and Wa arc marginalized over all other parameters. The solid 
line is for 5-band photometric redshift survey, the dashed line for 
a 9-band survey and the dot-dashed line for 17-band survey. Note 
we assume shapes are not measurable beyond Zmax = 1.5. 

6 SYNERGY OF DARK ENERGY 
EXPERIMENTS 

Here we present the results of comparing and combining 

3D weak lensing with other dark energy probes. Combining 
probes, which use different cosmological effects to measure 
dark energy either through the growth of structure or ge- 
ometric effects, will allow for cross checks to made. These 
cross checks may illuminate possible discrepancies between 
the two effects which could be important in determining 
the nature of dark energy. 3D weak lensing probes both 
the growth of structure via the matter power spectrum, and 
geometry through the lensing effect and the matter power 
spectrum. 

6.1 Comparing and combining with CMB, BAO 
and SNIA experiments 

We consider the CMB, SNIa and BAO experiments as de- 
scribed in Section 5.1. In Figures 12, 13 and 14 the dark thin 

ellipse is the CMB constraint; the small lightest gray ellipse 
is the lensing constraint; the darker gray, almost vertical, 
broad ellipse is the BAO constraint; the very broad lighter 
gray ellipse is the SNIa constraint. 

Figure 12 shows the combined two parameter 1-cr 
(68.3%) contours for all the possible pairs of experiments 
considered. In comparison with the other methods the 3D 
weak lensing constraint provides the smallest independent 
constraint. In combination the marginal errors so not vary 
largely between the different pairs. The BAO and CMB pair 
combination provides the smallest marginal errors through 
the unique degeneracy of the BAO ellipse providing a small 
intersection with the CMB. The SNIa constraint alone is 
poor although in combination with the other dark energy 
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-2-10 1 2 

Figure 12. The combined marginal wo, Wa constraints for four 
individual experiments combined in pairs. The experiments are a 
darkCAM lensing experiment and a CMB 14-month Planck ex- 
periment, a BAO WFMOS experiment and a SNIa SNAP exper- 
iment. The dark ellipses in the diagonal panels is the CMB con- 
straint; the small, light gray, ellipses along the top row of panels 
are the lensing constraint; the second lightest gray ellipses in the 
right-hand panels are the SNIa constraint; the broad darker gray 
ellipses in the top middle, right-hand middle and middle diagonal 
are the BAO constraint. The small white central ellipses are the 
combined constraints. 
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Figure 13. The combined marginal wq, Wa constraints for four 
experiments, combined three at a time. The experiments are a 
darkCAM lensing experiment, a CMB 14-month Planck experi- 
ment, a BAO WFMOS experiment and a SNIa SNAP experiment. 




-2 2 -2 2 

w w 

a a 

Figure 14. The left hand panel shows the combined marginal wq, 

Wa constraints for all four experiments combined allowing for fully 
open models, the right hand panel shows the combined marginal 
«)(), Wa constraints for all four experiments combined, with the 
condition Qm+^v = 1 enforced. The experiments are a darkCAM 
lensing experiment and a CMB 14-montli Planck experiment, a 
BAO WFMOS experiment and a SNIa SNAP experiment. 



probes does provide an improvement on the marginal errors 
through the intersection of the constraints. The SNIa is also 
a purely geometric test, so that the combination with 3D 
weak lensing would provide an important cross check. 

We combine combinations of three experiments in Fig- 
ure 13. It can be seen, by comparing with Figure 12, that 
in adding further information from another experiment im- 
proves all the combined marginal errors. The largest im- 
provements are gained by adding 3D weak lensing, or the 
CMB, to the pair combinations that do not include these 
probes. The smallest marginal errors are achieved by com- 
bining 3D weak lensing with the CMB and BAO. In com- 
bining 3D weak lensing with SNIa and BAO the dark energy 
constraints are comparable, or better than, each pair com- 
bination that includes the CMB constraint. 

Finally the combination of all four of the dark energy 
probes is shown in Figure 14. By adding 3D weak lensing to 
the three experiment combination of CMB, BAO and SNIa 
the marginalized constraints are improved by a factor of 2. 
The marginalized constraints in the combination of all four 



of the probes considered are 
Awo = 0.035 

Awa = 0.094, (43) 

with a pivot redshift error of 

Aw{z = 0.43) = 0.0147. (44) 

Reducing the maximum £ to 2000 increases these errors by 
approximately 0.01 to Awo = 0.045 and Awa = 0.105. 

6.2 Complementary figures of merit and pivot 

redshifts 

An illustrative way to present the information of pivot red- 
shifts and the figure of merit of a dark energy probe, or 
combination of different probes, is to show how the figure 
of merit and the pivot redshift compare. We show this in 
Figure 15 by plotting the figure of merit for all the possible 
combinations of experiments against the pivot redshift of the 
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Table 4. Expected marginal errors on cosmological parameters from 3D analysis of proposed weak lensing surveys. Note here 9 bands 
refers to 5 optical bands plus 4 infrared. 
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Figure 15. The figure of merit and pivot redshift for various ex- 
perimental combinations. The combinations are labeled as L=3D 
weak Lensing, B=BAO, S=SNIa, C=CMB. Combinations of let- 
ters represent combinations of experiments. 



combined constraint. In general the more experiments that 
are added in combination the smaller the figure of merit be- 
comes. As more experiments are added in combination the 
pivot redshift converges to a single value. In combination 
with other experiments the BAO constraint creates a high 



pivot redshift, this is due to the redshift of the nearest bin 
a,t z — 1. The 3D weak lensing constraint in combination 
creates a low z « 0.4 pivot redshift; this is due to the lens- 
ing pivot redshift dominating which is a symptom of the 
lensing signal maximizing at around that redshift. It can be 
seen from Figure 15 that there exists combinations, for ex- 
ample 3D weak lensing with the CMB (LC) and 3D weak 
lensing with BAO and SNIa (LSB), that have similar figures 
of merit but very different pivot redshifts. In using combina- 
tions such as these the w{z) evolution could be constrained 
to a high degree over a large redshift range. 



6.3 Future lensing surveys 

There are a number of current and planned imaging surveys 
for weak lensing which could be analyzed in 3D. The sur- 
veys vary in depth, areal coverage and number of bands, 
and illustrative errors are shown in the Table 4. Errors 
arc marginalized. The surveys considered are: the Canada 
France Hawaii Legacy Survey (CFHTLS; Semboloni et al., 
2006) which is ongoing; the VST (VLT Survey Telescope) 
public survey KIDS; SNAP (Supernova/Acceleration Probe; 
Aldering, 2005), and darkCAM on a 4 metre telescope. We 
show the errors achievable with darkCAM combined with 
various different experiments. BAO darkCAM refers to us- 
ing the photometric redshifts from darkCAM to measure 
BAO. VST-KIDS and CFHTLS have been combined with 
a 4 year WMAP prior as Planck will not be contemporary 
with these surveys. Here 9 bands refers to a 5 band optical 
survey with 4 infrared bands as discussed in Section 5.5. 



© 0000 RAS, MNRAS 000, 000-000 



16 A. F. Heavens, T. D. Kitching, A. N. Taylor 



6.4 The effect of changing the fiducial dark 
energy model 

Here we investigate the effect of the assumed fiducial dark 
energy cosmology on the marginal errors from our Fisher 
matrix calculations. The assumed cosmology has been a cos- 
mological constant model with wo = — 1 and Wa = 0. Figure 
16 shows the two parameter 1-cr (68.3%) contours for various 
dark energy models in the (wo, Wa) plane fully marginalized 
over other parameters. We consider two extreme examples, 
just allowable from current constraints: a SUGRA (Super 
Gravity) model proposed by Weller & Albrcclit (2002) rep- 
resented by Wo = —0.8 and Wa = +0.3; and a phantom 
model proposed by CaJdwell et al. (2003) with wo = —1.2 
and Wa = —0.3. Despite the marginal errors from the dark 
energy experiments alone changing, the combined marginal 
error on wq is largely unaffected by the assumed dark en- 
ergy model. The main difference is occurs on the error on Wa 
which increases for all methods as its value becomes more 
negative. This is duo to the fact that a negative Wa repre- 
sents a dark energy scenario in which the dark energy den- 
sity was less in the past (increasing into the future) ; so that 
the effect of dark energy on the expansion rate on observed 
galaxies (in the past) is less in these scenarios (and similarly 
the opposite effect for a positive Wa)- For the range of wq 
and Wa allowed by current constraints the marginalized er- 
rors presented here should be robust to the actual nature of 
dark energy. 
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Figure 16. The two parameter l-a (68.3%) contours for various 
assumed fiducial dark energy models, for a 10,000 square degree 
survey to a median depth = 0.7, with a 14- year Planck prior, a 
BAO WFMOS prior and a SNIa SNAP prior. The errors quoted 
are the one parameter l-a marginal errors on wq and Wa- The 
dark thin elhpse is the CMB constraint; the small hghtest gray 
ellipse is the lensing constraint; the darker gray, almost vertical, 
broad ellipse is the BAO constraint; the very broad light gray 
ellipse is the SNIa constraint. The small white central ellipses are 
the combined constraints. 



6.5 The effect of assuming flatness 

Here we present the effect of assuming flatness in the pa- 
rameter error estimation i.e. fi„ = 1 — f2,„. Figure 14 shows 
the two parameter l-a (68.3% contour) constraints for all 
four dark energy probes considered, see Section 5.1, in the 
(wo, Wa) plane with the assumption that the Universe is flat. 
The CMB 14-month Planck constraint does not considerably 
improve because the CMB puts a strong constraint on the 
overall geometry of the Universe, through the position of the 
first acoustic peak. The improvement in the SNIa constraint 
is most evident, since this dark energy probe is very sensi- 
tive to the overall geometry through the Hubble parameter 
(this is in agreement with Linder, 2005). The 3D weak lens- 
ing and BAO constraints also considerably improve, with 
the overall combined errors on the dark energy equation of 
state parameters being Awo = 0.029 and Awa = 0.089, 
a factor of 1.2 less than the constraints considering fully 
open models. In assuming flatness the different dark energy 
probes still have unique and complimentary degeneracies in 
the (wo, Wa) plane. The reduction in predicted errors, espe- 
cially in the SNIa, 3D weak lensing and BAO experiments, 
show that the assumption of flatness can have an affect on 
parameter error estimation (this is in agreement with the af- 
fect of this assumption on weak lensing tomographic meth- 
ods, see for example Knox, Song & Zlian, 2006). Given that 
some proposed dark energy models rely on modifications to 
the Friedmann equation in non-flat geometries it is prudent 
to calculate predicted parameter errors using fully open ge- 
ometries. 



7 CONCLUSIONS 

In this paper we have presented a 3D weak lensing spectral 
method suitable for high-£ studies, and investigated how 
well 3D weak lensing surveys could determine the equa- 
tion of state of dark energy. The accuracy which could 
be achieved if systematic errors can be controlled is im- 
pressively high, provided the surveys are analyzed in 3D: 
marginal statistical errors of Awq — 0.108, on the current 
value of w = p/(p(?), and its evolution Wa constrained to 
Awa = 0.397 are possible with a 10,000 square degree sur- 
vey in 5 bands to a median source depth of Zm = 0.7. At 
a pivot redshift oi z = 0.37 such an experiment could con- 
strain w{z) to Aw{z = 0.37) = 0.0175. Such a survey is 
possible with darkCAM, in conjunction with data from the 
Planck satellite. Even without Planck, the accuracy from 
3D weak lensing alone is still impressively high, and better 
than any other dark energy probe considered on its own. The 
fact that the physics of 3D weak lensing is well-understood, 
combined with the small statistical error forecasts, makes 
3D weak lensing a formidable prospect for advancing cos- 
mology in the next decade. The errors on w are comparable 
to, but a little better than, predictions from tomography (Hu 
and Jain, 2004; Ishak, 2005). The constraints on w{z) at the 
pivot redshift and the figure of merit, Aw{zpivot) * Awa, of 
the experiments considered were also discussed. 

We have investigated optimizing a wide field survey to 
measure the equation of state parameters wo and Wa and 
found an optimal survey strategy oi Zm = 1-0 covering 2400 
square degrees for a 5 optical band survey. We found that in- 
creasing the number of optical bands to 9 or 17 makes little 
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difference to the marginal errors when the 3D weak lensing 
result is combined with a Planck prior. The effect of includ- 
ing infrared bands in a wide field survey was investigated 
by varying the photometric redshift error, it was found that 
adding 4 infrared bands to a 5 band optical survey improves 
the marginal constraints on wo slightly from Awo = 0.108 
to Awo = 0.097. 

Three alternative dark energy probes were considered: a 
Planck CMB experiment; a WFMOS BAO experiment and 
a SNAP SNIa experiment. All possible combinations of ex- 
periments were considered and the figure of merit and pivot 
redshifts of the combinations shown. In such a competitive 
environment 3D weak lensing places strong constraints on 
the dark energy parameters and in combination with other 
experiments provides a unique degeneracy in the (u'o.uia) 
plane which is manifest as a strong constraint at a particu- 
lar pivot redshift. 

We have addressed the issues of biased photometric red- 
shift estimates (e.g. Ma, Hu and Huterer, 2005) and show 
that the method is relatively insensitive to this. Wo also 
investigated the effect that a sample of outliers, with poor 
photometric redshift estimates, would have on the predicted 
marginal errors. The effect of outliers on the marginal error 
of Wo is small although the way in which such a sample is 
treated is important. 

We have not considered errors due to the intrinsic align- 
ment of galaxies (Heavens, Refregior and Heymans, 2000; 
Croft and Mctzler, 2000; Catolan, Kamionkowski and Bland- 
ford, 2001; Crittenden et al., 2001; Jing, 2002), as these may 
be reduced to a negligible level by removing pairs which 
are close in photometric redshifts (Heymans and Heavens, 
2003; King and Schneider, 2002). This procedure has al- 
ready been demonstrated in the analysis of the COMBO-17 
data (Heymans et al., 2004). We have also not addressed 
other issues of systematics, such as optical distortions, or 
possible alignment of foreground galaxies with shear (Hirata 
and Seljak, 2004), which may be reduced using techniques 
such as template fitting (King, 2005). Nevertheless, the fact 
that the statistical errors are very small is very encourag- 
ing. Clearly to achieve the accuracies quoted here is going to 
be a formidable challenge for control of systematics, but at 
least the statistical error forecasts arc small enough that the 
promise of accurate measurement of the equation of state of 
dark energy may be realized. 
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